■  .... 


mm 


MC  FILE 


fib 

COP)  ^ 


OFFICE  OF  NAVAL  RESEARCH 


Contract  N00014-86-K-0043 


TECHNICAL  REPORT  No.  73 


Light  Scattering  by  a  Phase  Conjugator  in  the  Four- Wave 
Mixing  Configuration 


Henk  F.  Arnoldus  and  Thomas  F.  George 


Prepared  for  Publication 


Journal  of  Modern  Optics 


Departments  of  Chemistry  and  Physics 
State  University  of  New  York  at  Buffalo 
Buffalo,  New  York  14260 


June  1988 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the 
United  States  Government. 

This  document  has  been  approved  for  public  release  and  sale; 
its  distribution  is  unlimited. 


DTIC 

electei 

JUL  2  219881 


D 


t 


f 

It; 


V 

l 

I 


ft 


—  ■  ” k  VA\V.V.V_V\,V> 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


la.  REPORT  SECURITY  CLASSIFICATION 
Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 
UBUFFALO/DC/88/TR-73 


16.  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 
Approved  for  public  release;  distribution 

unlimited 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6b.  OFFICE  SYMBOL 
(If  applicable) 


8b.  OFFICE  SYMBOL 
(If  applicable) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 
Depts.  Chemistry  &  Physics 
State  University  of  New  York 


6c  ADORESS  (City,  State,  and  ZIP  Code) 

Fronczak  Hall,  Amherst  Campus 
Buffalo,  New  York  14260 


8a.  NAME  OF  FUNDING  /  SPONSORING 
ORGANIZATION 

Office  of  Naval  Research 


8c  ADORESS  (City,  State,  and  ZIP  Code) 

Chemistry  Program  program  project  itask 

800  N.  Quincy  Street  element  no.  no.  I  no 

Arlington,  Virginia  22217  _ I _ 

11.  TITLE  (Include  Security  Classification) 

Light  Scattering  by  a  Phase  Conjugator  in  the  Four-Wave  Mixing  Configuration 


12.  PERSONAL  AUTHOR(S) 


7a.  NAME  OF  MONITORING  ORGANIZATION 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 
Chemistry  Program 
800  N.  Quincy  Street 
Arlington,  Virginia  22217 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUM8ER 

Contract  N00014-86-K-0043 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROJECT 

TASK 

NO. 

NO 

WORK  UNIT 


Henk  F.  Arnoldus  and  Thomas  F.  George 


13a.  TYPE  OF  REPORT 


13b.  TIME  COVERED 
FROM  _____  TO 


14.  DATE  OF  REPORT  (Year,  Month.  Oay) 
June  1988 


16.  SUPPLEMENTARY  NOTATION 


Prepared  for  publication  in  Journal  of  Modern  Optics 


COSATI  CODES 


■>  .  '! 

(.  n  /»'  r ••  A- 


'h  f  A  - 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  ahd  identity  by  block  number) 
IGHT  SCATTERING;  -  > PHASE+MATCHING  RESONANCES 


PHASE  CONJUGATOR 


cessary  ahd  identity  by  block  number) 

PHASE+MATCHING  RESONANCES, 

FINITE  INCIDENCE  ANGLE 
SPECTROSCOPIC  APPLICATIONS 


FOURfrWAVE  MIXING  •  SPECTROSCOPIC  APPLICATIONS 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Reflection  of  travelling  and  evanescent  plane  waves  by  a  four-wave  mixing  phase 
conjugator  is  studied  in  detail.  No  restrictions  are  imposed  on  the  nonlinear  interaction 
strength,  the  angle  of  incidence  or  the  frequency  mismatch  between  the  pump  beams  and  the 
incoming  waves.  We— only  assume— that)  the  incident  field  is^weak  compared  to  the  pump  fields, 
which  justifies  a  classical-field  description  of  the  pumps.  The  wave  vectors,  amplitudes 
and  phases  for  the  various  waves  are  evaluated,  without  the  slowly-varying  amplitude 
approximation.  Familiar  phase-matching  resonances  for  certain  values  of  the  interaction 
length  are  recovered,  and  in  addition  strong  resonances  are  found  if  the  angle  of  incidence 
is  finite,  and  the  incident  light  is  not  in  perfect  resonance  with  the  pumps.  The  latter 
resonances  appear  at  the  transitions  from  a  travelling  to  an  evanescent  wave.  The 
significance  of  finite  angles  of  incidence  and  evanescent  waves  for  spectroscopic 
applications  is  pointed  out.  '  n 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
53  UNCLASSIFIEO/UNLIMITED  S3  SAME  AS  RPT 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 
Dr .  David  L.  Ne 


DD  Form  1473,  JUN  86 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

□  otic  USERS  Unclassified 


22b  TELEPHONE  (Include  Area  Code)  22 c.  OFFICE  SYMBOL 
(202}  696-4410 


Previous  editions  are  obsolete.  SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


'WWW  vvwwirv’iru  iru 


i  * 

Journal  of  Modern  Optica,  in  press 
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State  University  of  New  York  at  Buffalo 
Buffalo,  New  York  14260,  USA 

Abstract 

Reflection  of  travelling  and  evanescent  plane  waves  by  a  four -wave  mixing 
phase  conjugator  is  studied  in  detail.  No  restrictions  are  imposed  on  the 
nonlinear  interaction  strength,  the  angle  of  incidence  or  the  frequency 
mismatch  between  the  pump  beams  and  the  incoming  waves .  We  only  assume  that 
the  incident  field  is  weak  compared  to  the  pump  fields,  which  justifies  a 
classical-field  description  of  the  pumps.  The  wave  vectors,  amplitudes  and 
phases  for  the  various  waves  are  evaluated,  without  the  slowly- varying 
amplitude  approximation.  Familiar  phase -matching  resonances  for  certain  values 
of  the  interaction  length  are  recovered,  and  in  addition  strong  resonances  are 
found  if  the  angle  of  incidence  is  finite,  and  the  incident  light  is  not  in 
perfect  resonance  with  the  pumps .  The  latter  resonances  appear  at  the 
transitions  from  a  travelling  to  an  evanescent  wave.  The  significance  of 
finite  angles  of  incidence  and  evanescent  waves  for  spectroscopic  applications 
is  pointed  out. 


PACS:  03 . 50 . De ,  41.10.Hv,  42.65.Hw 


1 .  Introduction 


Among  the  many  ways  [1-3]  of  generating  a  phase -conjugated  signal  with 

respect  to  a  reference  signal,  the  technique  of  four-wave  mixing  is  the  most 

promising  from  an  experimental  point  of  view  [4-9],  A  nonlinear  crystal  (like 

BaTiO^)  or  liquid  (typically  CS2)  is  irradiated  by  two  counterpropagating 

strong  laser  beams  (the  pumps)  with  intensity  I.  A  third  incident  (weak)  field 

then  couples  to  the  pump  fields  through  the  third-order  susceptibility 

and  the  result  is  an  electric  polarization  of  the  medium,  which  is  proportional 
(3) 

to  x  I  and  the  electric  field  component  of  the  weak  field.  This  induced 
polarization  then  emits  radiation  which  propagates  out  of  the  crystal.  Under 
certain  conditions  this  generated  fourth  wave  is  the  phaase-conjugated,  or 
time -reversed,  replica  of  the  incident  field.  Production  of  phase -conjugated 
radiation  is  of  great  practical  importance  in  optical  engineering,  because  it 
provides  a  method  for  correction  of  wavefront  distortions. 

In  most  applications  the  weak  field  is  a  nearly-monochromatic  plane  wave 

with  well-defined  polarization,  the  angle  of  incidence  on  the  crystal  is  almost 

(3) 

zero  (usually  a  few  degrees),  and  the  coupling  constant  7  a.  x  I  is  very 
small.  For  this  configuration  the  generation  of  phase -conjugated  waves  is 
well -understood  [10-19].  There  are,  however,  conceivable  applications  in  which 
these  conditions  do  not  hold.  It  has  been  predicted,  for  instance,  that  the 
lifetime  of  an  atom  in  the  neighborhood  of  a  phase  conjugator  (PC)  is  infinite, 
as  a  result  of  the  fact  that  an  ideal  PC  focusses  the  emitted  fluorescence 
exactly  back  on  the  atom  [20,21],  Consequently,  the  spectroscopic  linewidth  of 
the  atomic  transition  under  consideration  would  be  zero,  which  can  have  a  great 
impact  on  frequency  standards.  These  predictions  were  derived  under  the 
assumption  of  perfect  phase  conjugation.  Since  phase  conjugation  is  equivalent 
to  time  reversal,  it  is  obvious  that  ideal  PC's  cannot  exist  (violates 
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causality).  Nevertheless  it  can  be  anticipated  that  realistic  PC's  can 
possibly  be  utilized  to  manipulate  linewidths  over  a  large  range  (in  contrast 
to  the  situation  of  atoms  near  a  metal  surface,  where  the  change  in  lifetime  is 
at  best  a  factor  of  two) . 

Emitted  dipole  radiation  (fluorescence)  by  an  atom  in  the  vicinity  of  a  PC 
has  plane-wave  components,  which  are  incident  on  the  surface  under  every  angle 
of  incidence.  Besides  that,  a  dipole  field  has  evanescent  components 
(exponentially-decaying  waves)  [22],  and  the  radiation  is  not  monochromatic. 
Furthermore,  for  contemporary  high-power  lasers  the  interaction  parameter  7  a  I 
is  not  necessarily  small.  In  this  paper  we  present  a  general  treatment  of  the 
scattering  of  travelling  and  evanescent  waves  by  a  four-wave  mixing  PC,  without 
restrictions  on  the  interaction  strength,  frequency,  polarization  or  angle  of 
incidence. 


A  nonlinear  transparent  crystal  (x^  **  0,  x^*^  "  °)  occupies  the  region 

0>z>  -  A,  A>0,  in  an  xyz  Cartesian  coordinate  frame,  and  the  regions  z  >  0 

and  z  <  -A  are  empty  space.  Two  counter -propagating  pump  beams  with  intensity 

I  and  frequency  u  (called  the  setting  frequency  of  the  PC)  illuminate  the 

(3) 

medium.  Then  the  complex-valued  coupling  parameter  is  given  by  7  oc  x  1-  We 

(3) 

shall  assume  that  I  is  constant  (no  depletion  of  pumps)  and  that  x  is 
frequency  independent.  This  means  that  the  incident  probe  field  must  have  a 


bandwidth  around  w  which  is  smaller  than  the  frequency  width  of  x' 


In  this 


fashion  we  can  avoid  complicated  notations,  but  it  is  straightforward  to  retain 

the  frequency  dependence  of  7  if  necessary  [23].  Furthermore,  we  assume  that 

(3) 

the  tensorial  nature  of  x  is  irrelevant,  which  can  always  be  managed  by  a 


proper  choice  of  geometry. 
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The  radiation  field  shall  be  represented  by  its  electric  and  magnetic 
components  E(r,t)  and  B(r,t),  respectively,  and  charges  and  currents  by  a 
polarization  density  P(r,t).  It  is  advantageous  to  adopt  a  Fourier  transform 
with  respect  to  time, 


dt  e  E(r,t)  ,  w  real, 


(2.1) 


and  since  E(r,t)  is  real  we  have 


£(r ,  -w)  -  £(r,<j)* 


(2.2) 


Similar  notations  apply  to  B(r,t)  and  P(r,t).  In  the  Fourier  domain  Maxwell's 
equations  read 


V-(£  +  p/f  )  -  0  ,  V-fi  -  0  , 


(2.3) 


V  x  £  —  iwfi 


,  V  x  6  -  ^  (£  +  P/tQ)  , 


which  should  hold  for  all  r  and  w. 


In  the  regions  z  >  0  and  z  <  -A  the  polarization  density  P  is  zero. 

Inside  the  nonlinear  medium  the  f’(r,w)  is  proportional  to  the  electric  field  at 
a  different  frequency.  Explicitly  [24],  this  has  the  fora 


TWSPnnmrH  THmuntx  u<  AT  w»n 


P(r,w) 


«q7*  E(r,w-2o>)  ,  w  >  0 


«q7  £(r,2aH-w)  ,  u>  <  0 


(2.4) 


where  the  field  £  does  not  include  the  two  pump  fields.  These  are 
parametrically  accounted  for  by  7  a  I.  On  the  surfaces  z  -  0  and  z  -  -A, 
Maxwell's  equations  (2.3)  imply  the  usual  boundary  conditions. 


3 .  Dispersion  relation 

Before  we  can  solve  the  scattering  problem  for  travelling  and  evanescent 
waves  by  this  PC,  we  have  to  establish  the  fundamental  plane -wave  solutions 
which  are  supported  by  the  medium.  To  this  end  we  first  notice  that  the 
polarization  density  ?(r,u)  from  Eq.  (2.4)  couples  positive  and  negative 
frequencies.  If  we  denote  by  »  w  >  0  a  fixed  positive  frequency,  then 
P(£  ,u>^)  is  determined  by  the  electric  field  component  with  frequency 

u>2  ”  "  2w  (3.1) 

The  polarization  at  this  negative  frequency  =  -w  is  then  proportional  to  the 
electric  field  at  2w  +  according  to  Eq.  (2.4).  Hence  the  nonlinear 

interaction  couples  positive  and  negative  frequencies  and  in  pairs. 

Consequently , for  a  fixed  Maxwell's  equations  (2.3)  constitute  essentially  a 
set  of  eight  equations,  which  have  to  be  solved  simultaneously. 

The  third  Maxwell  equation  can  be  written  as 


8(r,w)  -  x  £(r,«)  , 


(3.2) 


* 


S3 


and  therefore  8  is  known  as  soon  as  we  have  found  £ ,  both  for  and  .  Then, 

V*fi  -  0  is  automatically  satisfied.  Furthermore,  we  notice  that  £  is 
proportional  to  £,  although  at  a  different  frequency,  and  so  the  first  Maxwell 
equation  is  certainly  obeyed  if 


V«£(r,u>)  -  0  , 


(3.3) 


for  every  w.  Next  we  substitute  Eqs.  (2.4)  and  (3.2)  into  the  fourth  Maxwell 
equation,  which  yields  the  set  of  coupled- wave  equations 


9  19 

{?  +  (-£)z)  £(r.Wl) 


9  9  9 

(V  +  (-^)Z)  fi(r,«2) 


-t*(-7>2  £(£,«2)  . 


y  9 

■7  (.-%)  £<£.«!)  , 


(3.4) 


(3.5) 
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: 
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W. 
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for  the  electric  field.  Equations  (3. 3) -(3. 5)  are  the  basic  relations  for  a 
PC.  The  first  one  states  that  the  fields  are  transverse,  and  (3.4)  and  (3.5) 
show  that  a  wave  with  a  positivie  frequency  couples  to  a  wave  with  a 
negative  frequency  and  vice  versa.  Therefore,  a  positive-frequency  field 
generates  a  negative- frequency  field,  which  is  the  essence  of  a  phase 
conjugator . 

As  a  plane -wave  solution  we  try 


-  h  e 


ik  *r 


5 

I, 

I 


(3.6) 


ik  •  r 

^a(J’w2)  “  ^a~a  e  ~ 


(3.7) 


s 

» 

$ 
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Then  Eqs .  (3.3)  -  (3.5)  give 


k  *E  -  0 
-a  -a 


-  k2(l  +  7*^)  -  kV(l  +  7/ua) 


where  we  introduced  the  wave  number 


(3.8) 


(3.9) 


k  —  w^/c 


and  the  dimensionless  detuning  parameter 


p  -  2w/w^  -  1 


(3.10) 


Equation  (3.3)  states  that  plane  waves  in  a  PC  are  transverse,  and  the  last 
equality  in  Eq.  (3.9)  gives  an  equation  for  the  amplitude  ratio  rj  between  the 

c l 

and  the  component.  Because  Eq.  (3.9)  is  quadratic  in  rj^,  it  admits  two 
solutions.  For  reasons  which  will  become  clear  in  due  course,  we  choose  the 
solution 


„a  -  {p2-l-S((p2-l)2  +  47oP2)W)/27*  , 


(3.11) 


7  -  yne  v  ,  yn  >  0  ,  <A  real 


(3.12) 


(3.13) 


S  -  sgn(w-ui1  )  -  sgn(p-l) 


The  solution  corresponding  to  the  second  root  for  will  be  written  as 


?b(H*wl>  -  rth ^  e 


lb(H-w2)  “  h  e 


*VE 


which  obeys  Maxwell's  equations  if 


V*b  -  0 


k£  -  k2p2(l  +  7^b)  -  k2(l  +  7*/nb) 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


Of  the  two  possible  solutions  for  we  have  to  take 


r/b  -  (1  -  P2  +  SUP2-!)2  +  47qP2)H)/27P2 


(3.18) 


With  this  convention  we  have  n  -*0,  rj,  -»  0  if  7  -*  0,  and  the  a  and  b  solutions 
become  (uncoupled)  and  waves,  respectively,  in  this  limit.  For  7^0  the 
parameters  17 ^  and  determine  the  relative  strengths  of  the  coupled  waves  with 
the  complementary  frequency,  which  are  excited  by  the  four-wave  mixing  process. 
Furthermore,  we  notice  that  the  coupled  waves  £  (r,w, )  and  £  (r,u>„)  have  the 
same  wave  vector  k^,  which  implies  a  perfect  phase  matching  between  these  two 
waves  at  different  frequencies.  The  same  holds  for  the  b-solution. 


9 


»«*  *•#:*«!**** 


Now  we  can  substitute  the  expressions  for  i)  and  17,  into  Eqs .  (3.9)  and 

a  d 

(3.17),  which  gives  the  wave  numbers  and  (up  to  a  minus  sign).  Notations 
can  be  simplified  considerably  by  the  introduction  of  the  quantity 


£  -  H(p2+l-6((p2-l)2  +  47q pV) 


Then  we  find  in  terms  of  e 


(3.19) 


k  £  ,  k^  -  kz(p  +1-0 


(£-l)/7*  ,  r?,  -  (l-£)/7p 


(3.20) 


(3.21) 


2  2 

If  the  medium  would  be  an  ordinary  dielectric,  we  would  also  have  k^  -  k  e  ,  but 
with  e  as  the  dielectric  constant.  An  important  difference  is  that  in  Eq. 
(3.19)  the  £  depends  explicitly  on  the  frequency  through  the  parameter  p. 
This  frequency  dependence  is  a  genuine  geometrical  effect,  as  it  follows  from 
the  mechanism  of  four-wave  mixing  (rather  than  from  a  frequency  dependence  of 
\  which  we  have  suppressed).  Therefore,  Eq.  (3.20)  gives  the  fundamental 
relations  for  the  two  branches  of  the  dispersion  curve  for  a  four -wave  mixing 
PC.  This  universal  dispersion  relation  is  plotted  in  Fig.  1.  We  remark  that 
£ ,  as  defined  in  Eq.  (3.19),  is  real  (possibly  negative).  If  the  frequency 

(3) 

dependence  of  x  would  have  been  retained,  then  £  could  also  have  an 
imaginary  part.  Furthermore,  we  notice  that  e  is  discontinuous  across  the 
resonance  w,  -  w,  or  p  -  1,  due  to  the  appearance  of  5. 


Incident  field 


4. 

A  given  external  field  illuminates  the  surface  z  -  0  of  the  PC.  Since 
almost  every  field  can  be  expanded  in  plane  waves,  and  since  Maxwell’s 
equations  for  this  PC  are  linear  (even  though  the  four-wave  mixing  process  is 
not) ,  it  suffices  to  consider  an  incident  field  of  the  form 


£.  (r ,w1 ) 
-inc  -  1 


E. 

-me 


ik*  r 
e 


(4.1) 


defined  in  the  region  z  >  0.  The  corresponding  6  field  follows  from  Eq.  (3.2) 
In  Eq.  (4.1)  the  polarization  and  amplitude  E^nc>  and  the  wave  vector  k  are 
arbitrary,  with  the  only  restrictions 


k*E. 

-  -inc 


k*k  -  (Wl/c)‘ 


(4.2) 


according  to  Maxwell's  equations  in  z  >  0. 

With  the  unit  vector  e z  as  the  normal  to  the  surface,  we  can  decompose  k 
into  its  parallel  and  perpendicular  components  with  respect  to  the  xy-plane. 
We  write 


k  -  kii  +  k  e 

-II  z~z 


(4.3) 


and  similarly  for  any  other  vector  quantity.  Combining  this  with  Eq.  (4.2) 
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The  quantity  kZ  is  a  given  positive  number  for  a  fixed  ,  but  the  components 
k. ||  and  kzez  of  k  can  be  anything,  as  long  as  restriction  (4.4)  is  satisfied. 
For  most  practical  cases  it  is  sufficient  to  consider  only  real-valued  vectors 
kii ,  and  this  will  be  assumed  from  now  on.  We  shall  regard  the  quantities  k  - 


u>^/c  >  0  and  k|  as  given,  but  arbitrary.  Then,  the  right-hand  side  of  Eq. 

(4.4)  is  a  given  real  number,  and  there  are  two  possible  solutions  for  kz- 
Because  the  external  field  is  generated  by  sources  in  the  region  z  >  0,  we  have 


to  choose  the  causal  solution,  which  is 


7k2 -kjf  ,  if  k  >  kit 


■  i-/k jj - k2  ,  if  k  <  k|.  , 


(4.5) 


where  ky  -  /k|| *k||  >  0.  For  k  >  k|  the  kz  is  negative  and  real,  corresponding 
to  an  incident  travelling  plane  wave  from  the  region  z  >  0.  In  the  case  k  < 


kjj ,  kz  is  imaginary,  and  the  root  is  chosen  in  such  a  way  that  the  wave  decays 
exponentially  to  zero  in  amplitude  in  the  negative  z-direction.  This 
evanescent  wave  decays  in  the  direction  perpendicular  to  the  surface,  and 
travels  along  the  surface  in  the  k« -direction. 


5.  Fields 


Hie  incident  field  £^nc(r,w^)  induces  a  nonlinear  polarization  in  the 
medium,  which  in  turn  emits  radiation  according  to  the  coupled-wave  equations 
(3.4)  and  (3.5).  This  radiation  travels  out  of  the  crystal  and  gives  rise  to 


reflected  and  transmitted  waves  by  the  layer.  It  will  turn  out  that  the  fields 


.1 


everywhere  in  space  can  be  expressed  in  plane  travelling  or  evanescent  waves, 
depending  on  k  and  k| ,  Every  wave  a  will  therefore  have  a  spatial  dependence 
of  E^exp ( ik^ • r ) .  Transversality  then  requires 


k  *E  -  0  , 

~a  -a 


(5.1) 


for  every  wave.  Since  the  waves  must  be  matched  across  the  planes  z  -  0  and  z 
-  -A  with  the  aid  of  boundary  conditions,  which  should  hold  for  all  r  in  z  -  0 
and  z  -  -A,  all  wave  vectors  must  have  the  same  parallel  component.  Therefore, 
we  have 


k  -  kit  +  k  e 
-a  ~||  a,z-z 


(5.2) 


and  only  k  remains  to  be  determined.  On  the  other  hand,  the  dispersion 
Ot  j  z 

2  2  2 

relations  in  vacuum  and  in  the  PC  fix  the  value  of  k  -  kii  +  k 

a  II  a ,  Z 

Consequently,  the  only  freedom  we  have  left  is  the  choice  of  the  sign  of  ^ . 
In  the  regions  z  >  0  and  z  <  -A  this  sign  is  determined  by  the  requiremen  that 
the  waves  must  emanate  from  the  PC,  in  the  same  way  as  we  found  the  sign  of  k^ . 
For  the  fields  inside  the  PC  there  is  no  a  priori  way  to  fix  the  signs  of  the 
z -components  of  the  wave  vectors,  and  therefore  we  have  to  retain  all  possible 
combinations.  We  shall  only  write  down  the  expressions  for  the  electric 
fields.  Then,  the  magnetic  fields  can  be  found  from  Eq.  (3.2). 


5 . 1  Region  z  >  0 

From  the  arguments  above  it  follows  that  the  most  general  pi 
solution  in  the  region  z  >  0  is  given  by 


ane-wave 


% 


VA’/.AA'VA',v;..vjv;.\va 
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>  4 -4 


£(r,w, )  -  E.  e  +  E  e 
-  -  1  -i.nc  -r 


ik  •  r 

6(r,«,)  -  E _  e  ~pc  ~ 


ik  *r 


*■*  v  *•  t  _ 

-  —  Z  -pc 


(5.3) 


(5.4) 


At  frequency  there  is  only  one  other  possible  wave,  which  is  the  specularly  - 


reflected  r-wave  with 


k  -  k  ,  k  -  -k 
r  r ,  z  z 


(5.5) 


Although  this  field  resembles  the  reflected  wave  by  an  ordinary  dielectric,  it 
is  here  entirely  generated  by  the  four-wave  mixing  process  (we  have  set  - 
0) .  At  o>2  we  have  the  phase -conjugated  pc-wave  with 


kpc  “  "VC  >  °  1 


,2  J-  ,2  ,22 

k  -  k  -  kii  -  k  p 

pc ,  z  pc  (I 


(5.6) 


(5.7) 


Because  w_  is  negative,  the  pc-wave  travels  in  the  -k  direction  if  k  is 
2  r  -pc  pc , z 

real.  In  the  case  of  an  evanescent  pc-wave  the  wave  should  die  out  in  the 
positive  z-direction.  Consequently,  the  root  should  be  taken  as 


:  >  z 


-ykV-kn  ,  if  kp  >  kii 


i7k?-k2p2  ,  if  kp  <  kii 


(5.8) 
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Then  it  remains  to  determine  E 


and  E 
-pc 


5 . 2  Region  0  >  z  >  -A 

Inside  the  PC  the  fields  are  combinations  of  a-  and  b-solutions,  and  we 


have 


x  ik  *r  ik  *r  ik?”*r  ikT*r 

£(r,w1)  -  r  e  ~a  "  +  Ea  e  ~a  “  +  e  ~  +  e  "  ,  (5.9) 

ik+*r  ik  «r  ilc"*r  ikl*r 

£(r,w»)  -  r/  E+  e  ~a  ~  +  n  E'  e~a“  +  E^e  b~  +  E‘e  (5.10) 

•  ~  Z  a**S  ~  D  ~D 

2  2 

The  values  of  k&  and  k^  are  given  in  Eq.  (3.20),  and  for  the  z-components  of 
the  wave  vectors  we  write 


+ 

k~ 


a.z 


-  ±  k. 


kb, 


z 


±  k. 


with 


k 


2 

1 


k 


2 

2 


(5.11) 


(5.12) 


(5.13) 


(5.14) 


The  roots  are  taken  (arbitrarily)  as 


-A\-V\  ,  if  ^  >  k|  , 


kl  “ 


(5.15) 


-^k]rka  ’  if  ka  <  k5  > 


-y^-kf  ,  if  k*  >  k| 


(5.16) 


W_kb  ,  if  k£  <  kf 


5 . 3  Region  z  <  -A 

The  only  waves  which  travel  or  die  out  in  the  negative  z -direction  are 


!(r,u>i)  -  Ec  e 


ik*r 


(5.17) 


"  ?ni  e 


ikni  *  - 


(5.18) 


a>2 ,  which  has 


Vc  ™  -  Vc 
ni , z  pc , z 


(5.19) 
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where  the  wave  vector  of  the  transmitted  (t)  wave  is  the  same  as  the  incident 
wave  vector.  Furthermore,  there  is  possibly  a  nonlinear  (ni)  wave  of  frequency 


n 


According  Co  Eq .  (5.2)  every  (complex-valued)  wave  vector  kQ  lies  in  the 
plane  of  incidence,  spanned  by  the  (real-valued)  vectors  ki|  and  e^.  The 
ampliude-polarization  vectors  Eq  can  be  decomposed  into  a  surface  (s)  polarized 
and  a  plane  (p)  polarized  component,  which  are  perpendicular  to  and  lie  in  the 
plane  of  incidence,  respectively.  Since  is  restricted  by  k^-E^  -  0,  the 
only  ambiguity  in  a  decomposition  along  unit  s-  and  p-polarization  vectors  is 
the  choice  of  the  phase  of  the  unit  vectors.  We  take 


-as  ”  kii 


x  £z} 


(6.1) 


e 

~ap 


)  , 


(6.2) 


and  it  is  easy  to  check  that  e  ,  e  and  k  /k  constitute  an  orthornormal  set 

-as  -ap  ~a'  a 

of  unit  vectors  for  every  a.  Notice  that  e^s  is  defined  as  independent  of  a 

and  is  real.  The  p-polarization  vector  can  be  complex.  For  all  waves  only  the 

2 

value  of  k^  (real)  is  prescribed  by  the  dispersion  relation,  and  we  take  the 
square  root  as 


A2  ,  if  k2  >  0 

cl  a 

r 

k  -  * 
a 

i7-k2  ,  if  k2  <  0 

cl  a 


(6.3) 


It  can  be  proven  that  every  wave  is  an  s(p)-wave  if  the  incident  wave  is  an 
s(p)-wave.  Therefore,  we  can  distinguish  between  two  cases  and  write 


for  every  a.  The  simplification  lies  in  the  following  relations  for  the 


i 

j  2 -components  of  E^: 


I 

E  -  0  ,  s  waves 

I  a, z 


E 

a ,  z 


p  waves 


(6.5) 


Then  the  Fresnel  coefficients  X  (Y  )  for  s(p) -waves  are  defined  as  the  ratio  of 

a  a  r 

E^  and  E^nc.  or  equivalently  the  amplitudes  E^  are  written  as 


E  -  X  E.  e 
—a  a  me  -as 


for  s -waves 


E  -  Y  E.  e 
-a  a  me  ~ap 


for  p -waves 


(6.6) 


The  amplitude  Einc  of  the  incident  field  is  a  given  quantity,  and  the  unit 
polarization  vectors  are  geometrically  determined.  Therefore,  knowledge  of 
X^  and  Y^  for  every  wave  a  determines  the  scattering  of  any  wave  by  the  PC. 

7 .  Solution 

Maxwell's  equations  (2.3)  state  that  at  the  boundaries  z  -  0  and  z  -  -A 

the  tangential  component  of  E,  the  normal  component  of  E  +  P/cq,  and  the  B 

field  must  be  continuous,  both  for  and  In  matching  the  fields  in  the 

three  regions  across  the  boundaries ,  we  can  evaluate  all  E  ' s ,  and  the  results 

-a 


can  be  expressed  in  terms  of  the  Fresnel  coefficients  and  for  s-  and  p- 
waves ,  respectively.  It  is  convenient  to  express  the  Fresnel  coefficients  in 
dimensionless  quantities,  rather  than  in  wave  numbers.  Besides  p  and  «  we 
introduce 


u  -  - kz/k  , 


(7.1) 


which  allows  us  to  write  for  the  parallel  components 


2  2  2 
k f,  -  k  (1-u  ) 


(7.2) 


For  a  travelling  incident  wave  u  is  real,  restricted  by  0  <  u  <  1,  and  u  equals 
the  cosine  of  the  angle  of  incidence.  If  the  incident  wave  is  evanescent,  then 
u  is  positive  imaginary.  Furthermore,  we  define  dimensionless  wave  numbers  by 


m  -  k  /k  ,  a  -  a ,  b ,  1 , 2 
a  cr 


(7.3) 


m  -  k  /k  , 
P  pc,z' 


(7.4) 


and  the  layer  thickness  d  in  units  of  a  wavelength  of  the  incident  radiation 


d  -  kA/2 n 


(7.5) 


Due  to  boundary  conditions  at  z  -  -A,  phase  factors  appear  which  can  be 
expressed  as 


.L  —  ^ - l_ 

*1  “ 


.L  .  O  — 


-2 


L<»n1 


<• .  .m  >J. M. /a*. ».t  •.|».|  >ai  f1|  «(i  M  «.i  ^r4.i  i<i  v«(  m  > 


(7.6) 


0,  -  2xdu  0,  -  2xdm 
3  4  p 


For  travelling  waves  these  phases  are  real,  and  for  evanescent  waves  they  are 
positive  imaginary. 

After  laborous  computations,  it  follows  that  the  Fresnel  coefficients  can 

be  expressed  in  terms  of  eight  dimensionless  parameters  x^,  y^,  i  =*  1 . 4  as 

follows : 


X+  -  x. 
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a  2 
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(7.7) 


Here,  the  parameters  and  y^  for  s-  and  p-waves,  respectively,  are  solutions 


of  the  linear  sets 


(7.8) 


where  the  matrices  P  and  Q  are  given  by 


fJa(mp'ml) 
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-fla(ml+m  p)e 


,'a(Vml) 

(m1-u)e 
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P  2 
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(7.9) 
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yi^e  (m1(m1+mp)  -  (1-u  ) ( (ma/p)  -1) ) 


m,  (m.+u)  -  (1-u  )  (m  -1) 
i  l  a 

,?a(ml(ml+IV  '  ((ma/p)2-l)  ) 

- 

e  1(m1(m1-u)  -  (1-u2) (m2-l) ) 

"  2  2 
r7ae  {m1(m1-mp)  -  (1-u  )((m J p)  -1) 


nb(m2(m2-u)  -  (1-u2) (m£-l) } 
m2(m2-mp)  -  (1-u2) ( (n^/p)2-!) 

■  i-^2  2  2 

r?be  {m2(m2+u)  -  (1-u  )(m£-l)} 

"  ^2  2  2 
e  (m2(m2+m  ) - (1-u  ) ((n^/p)  -1) ) 
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(7.10) 


The  sets  of  equations  from  Eq.  (7.8)  can  readily  be  solved  analytically,  but 
the  resulting  expressions  are  lengthy  and,  in  turn,  not  transparent. 
Numerically,  one  solves  the  sets  directly,  rather  than  inverting  the  matrices  P 
and  Q. 


Although  the  results  of  the  previous  section  apply  to  any  situation,  in 
many  practical  cases  the  solution  can  be  simplified  considerably  because  of 
restrictions  on  the  order  of  magnitude  of  various  parameters.  It  is 
elucidating  to  work  out  a  special  limit  in  order  to  reveal  the  fundamental 
structure  of  the  Fresnel  coefficients.  If  the  incident  field  is  a  visible 
narrow-bandwidth  laser  and  exactly  on  resonance  with  w,  then  the  relative 
detuning  is  of  the  order  of  |p-l|  —  10  Furthermore,  the  nonlinear  coupling 


parameter  7^  has  an  order  of  magnitude  of  10  *  10  even  for  very  strong  CW 

pump  fields.  In  this  section  we  consider  the  limit  7^  -*■  0  (weak  interaction) 
and  p  -  1  (resonance),  which  implies  7q  «  |u  |  and  | p - 1 |  «  |u  |.  First  we 
expand  the  matrices  P  and  Q  up  to  leading  order  in  p  -  1  and  7q,  and  in  the  end 
we  take  the  limits.  For  this  situation  the  amplitude  factors  are  related  by 


(8.1) 


and  q equals 


*  fj  1  ♦  cf  )2I 

70  70 


(8.2) 


which  is  not  necessarily  a  small  parameter.  Equation  (8.1)  expresses  that  the 

coupling  strength  between  the  two  fields  of  the  a- solution  equals  the  coupling 

strength  between  the  two  components  of  the  b-solution.  This  must  be  so  in  this 

2  2 

limit,  since  p  -  1,  7  -*  0  implies  k^  -  k^  and  consequently  the  two  branches  of 
the  dispersion  relation  coincide. 


8 . 1  Travelling  waves 

For  0  <  u  <  1  the  incident  field  is  a  travelling  wave.  In  lowest  order  we 
find  for  the  relative  wave  numbers 


m  ”  m.  —  m0  **  -u  , 
P  1  L 


ma  “  ”b  "  1  ’ 


(8.3) 


and  the  Fresnel  coefficients  are  found  to  be 


X*  -  Y*  -  XT  -  Y‘  -  X  -  Y  -  X  4  -  Y  .  -  0  , 

a  a  d  b  r  r  ni  ni 
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(8.5) 


This  limit  has  several  remarkable  features.  First,  the  Fresnel  coefficients 
for  s-  and  p-waves  are  identical,  and  therefore  the  scattering  process  is 
polarization  independent.  Second,  the  specularly- reflected  wave  disappears, 
and  hence  the  field  which  is  reflected  back  into  the  region  z  >  0  consists 
entirely  of  the  phase-conjugated  signal  with  respect  to  the  incident  beam 
(times  a  factor) .  Third,  the  structure  of  X^c  is  completely  determined  by 
phase  factors,  in  the  combination 


^  +  i>2  ~  -2ird(m1-m2) 


(8.6) 
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According  to  Eq .  (8.3),  we  have  -  a,  »  0  in  first  order  in  p-1  and  7q ,  which 
would  make  ^  +  xl> 0  and  thereby  Xpc  -  0.  However,  m^  -  m2  is  multiplied  by 
the  relative  layer  thickness  d.  For  an  interaction  region  of  a  few 
centimeters,  we  have  d  -  10^  and  +  ^2  *  even  in  lowest  order.  Since 
and  m2  are  the  relative  wave  numbers  (of  the  z-component)  of  the  a-  and  b-waves 
in  the  PC,  respectively,  we  conclude  that  the  phase -conjugated  signal  is 
brought  about  by  constructive  interference  between  the  a-  and  b -modes  of  the 


If  the  incident  field  is  in  very  close  resonance  with  the  setting  w  of  the 
PC,  then  we  have 


r)  -  -6e 
'a 


I ’a  I  "  1  • 


(8.7) 


(More  precisely:  if  { p - 1 {  «  .  Under  this  condition  the  X  reduces  to 


Xpc  “  '  ) ) 


and  the  PC  reflectivity  becomes  infinite  for 


(8.8) 
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^1  +  i>2  “  (2n+l)7r  ,  n  integer 


(8.9) 


which  is  the  famous  resonance  condition  [3]. 

Even  if  |  p - 1 1  is  not  much  smaller  than  7q,  the  denominator  of  X  ,  Eq. 
(8.5),  has  still  a  resonance  at  the  solution  of 
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which  leads  again  to  condition  (8.9).  With  6  as  the  angle  of  incidence,  we 
then  find  that  the  PC  is  resonant  for  layer  thicknesses 


d  -  (n+h) 


Ja-pY 
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y0 


n  -  0,1,2. . . 


(8.11) 


At  the  resonance  the  value  of  X  is  found  to  be 

pc 


2»7, 


pc 


i  -  b. 


1 2  1 


(8.12) 


and  in  between  resonances,  where  exp(i(i^+^2) )  “  1.  we  have  X^c  -  0.  This 
resonance  behavior  is  illustrated  in  Fig.  3. 

An  important  conclusion  is  that  when  a  PC  is  resonant  for  radiation  under 
normal  incidence  ( 6  -0),  it  is  off-resonant  for  radiation  which  illuminates 
the  surface  under  a  finite  angle.  In  spectroscopic  applications,  where  the 
incident  field  is  dipole  radiation,  all  plane-wave  components  strike  the  PC  at 
a  different  angle,  and  therefore,  this  device  cannot  operate  as  a  perfect  phase 
conjugator  for  the  entire  field. 
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8.2. 


Evanescent  waves 

For  an  evanescent  incident  wave  we  have 
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-m^  -  ~  u 


(8.13) 


and  the  Fresnel  coefficients  become 


and  the  same  expressions  hold  for  p-waves.  We  notice  that  there  is  no 

reflection  at  all  back  into  the  region  z  >  0.  Furthermore,  if  the  layer 

thickness  d  is  much  larger  than  the  penetration  depth  l/|u|  of  the  waves,  then 

the  fields  in  z  <  -A  also  vanish.  We  conclude  that  there  is  hardly  any 

reflection  of  evanescent  waves  in  the  limit  p  -*  1,  7^  -*  0.  In  Fig.  4  we 

compare  the  Fresnel  coefficients  Ix^)  for  travelling  and  evanescent  waves. 

For  evanescent  waves  the  nonlinear  interaction  region  is  limited  to  the 

penetration  depth,  which  is  a  few  optical  wavelengths.  The  PC  cannot  generate 

much  radiation  in  such  a  thin  layer,  which  explains  the  very  small  values  of  X 

and  X  in  this  case . 
pc 


9 .  Resonances 

In  the  previous  section  we  found  that  jx  j  acquires  extreme  values  if  the 
dimensionless  layer  thickness  d  is  related  to  the  angle  of  incidence  9,  the 


detuning  p  and  the  coupling  parameter  yn  according  to  Eq.  (8.11).  These 


resonances  appear  due  co  perfect  phase  matching  of  the  a-  and  b-waves  in  the 


PC,  as  expressed  by  Eq.  (8,10).  For  |p-l|  «  7q  we  have  |  rj  ^  |  -  1,  and  then 
|Xpc|  becomes  infinite,  as  follows  from  Eq.  (8.12).  Beside  these  interference 
resonances,  a  PC  has  a  different  kind  of  resonances  which  appear  if  we  allow 
the  angle  of  incidence  to  be  nonzero  and  the  waves  to  become  evanescent.  If  we 

solve  Eq.  (7.8)  for  the  eight  parameters  x^,...,x^,  y^ . y^,  then  the  general 

expression  for  every  parameter  is  a  3  x  3  determinant  (because  of  the  zero's  on 
the  right-hand  sides),  divided  by  det(P)  or  det(Q).  Resonances  then  occur  for 
values  of  p,  u,  j ^  and  d  at  which  det(P)  and  det(Q)  is  very  small.  For 
instance,  in  the  limit  of  section  8  we  have 


-  i  \f) 

det(P)  -  (2u)4(e  2  + 

for  travelling  waves,  and 

det(P)  -  (2u)4(1  +  |» 

Si 
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2) 2 


-i(01+02) 
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(9.1) 


(9.2) 


for  evanescent  waves.  Then  it  is  obvious  that  the  right-hand  side  of  Eq.  (9.1) 

has  a  minimum  if  the  phase-matching  condition  (8.10)  holds,  whereas  the  right- 

hand  side  of  Eq.  (9.2)  has  no  pronounced  minima. 

Without  proof  we  state  that  the  second  kind  of  resonances  can  appear  there 

where  one  of  the  generated  waves  turns  from  a  travelling  wave  into  an 

evanescent  wave,  or  equivalently,  at  the  branch  points  of  the  square  roots 

which  define  the  z-components  of  the  wave  vectors  (section  5).  From  the 

2 

expressions  of  section  5,  in  combination  with  relation  (7.2)  for  kj|  ,  we  find 
that  the  various  waves  are  evanescent  under  condition: 
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Then  the  resonances  are  located  atu  -0,u  -  1  -  p  ,  u  -  1  -  £  and  u  -  e  - 

2  2 

P  .  Most  obvious  is  the  case  u  -  0,  or  u  -  0,  for  which  det(P)  must  be  small 

according  to  Eqs.  (9.1)  and  (9.2).  We  notice  that  the  right-hand  sides  of 

(9.3)  depend  only  an  7 ^  and  p,  and  not  on  u.  For  fixed  7q  and  u,  the 

resonances  appear  at  a  certain  detuning  p  between  of  the  incident  field  and 

the  PC  setting  frequency  u>.  For  7q  and  p  fixed  we  can  regard  the  resonance 

conditions  as  an  equation  for  u  (angle  of  incidence  or  inverse  penetration 

depth)  at  which  |X  |  and  the  other  Fresnel  coefficients  have  sharp  peaks.  We 
pc 

notice  that  the  positions  of  the  resonances  are  independent  of  d,  in  contrast 
to  the  resonances  of  the  previous  section. 

Let  us  take  p  and  fixed,  and  consider  the  behavior  of  the  Fresnel 


coefficients  as  a  function  of  u.  Then  the  resonances  are  located  at  u  -  u 

where  the  u  's  are  solutions  of 
res 
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b- waves 

:  u 
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provided  that  the  equations  have  a  solution  in  the  range  of  u.  We  have 

suppressed  the  case  u  -  0,  since  the  corresponding  extrema  are  minima.  The 

2  2 

right-hand  sides  of  Eq.  (9.4)  are  real,  and  the  range  of  u  is  <  u  <1.  To 
see  the  physical  significance  of  the  resonance  conditions,  we  look  at 
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u  »  1  -  o  for  Che  pc-wave,  and  the  same  picture  will  hold  for  the  a-  and  b- 
res  r 

waves.  If  1  -  p  >1  there  is  no  solution,  and  the  Fresnel  coefficients  will 

2 

vary  smoothly  as  a  function  of  u.  For  0  <  1  -  p  <1  there  is  a  solution  with 

0  <  u2  <1  which  implies  0  <  u  <1  because  the  u-values  are  restricted  by 
res  ’  res 

0  <  u  <  1,  and  u  -  iv  with  v  >  0.  This  situation  corresponds  to  a  travelling 

incident  wave  with  u  -  u„  .  For  u  >  u  the  pc-wave  is  a  travelling  wave, 

res  res 

and  for  u  <  u  the  pc-wave  is  evanescent.  Exactly  on  the  transition  between 
res 

the  two  situations,  the  Fresnel  coefficients  have  a  sharp  resonance.  In  the 


case  that  1  -  /  <  0  we  write  urgs  -  i-vres.  and  the  solution  is  vrejj  -  ■ 

Then  the  incident  wave  is  evanescent,  and  the  pc-wave  is  evanescent  for  v  > 

v  and  travelline  for  v  <  v  We  notice  the  remarkable  fact  that  an 

res  °  res 

evanescent  incident  wave  can  be  reflected  by  the  pc  as  a  travelling  wave. 

In  Figs.  5  and  6  we  have  plotted  |Xp(J  for  a  travelling  and  an  evanescent 
incident  field,  respectively.  Curves  a  and  b  correspond  to  u^/u  -  1.05  and 
wx/w  -  0.95,  which  gives  p  -  0.905  and  p  -  1.105,  respectively.  We  chose  the 
value  0.05  for  parameter  7.  Then  the  solutions  of  Eq.  (9.4)  are  for  curves  a 
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and  for  curves  b 


pc-wave 
a -waves 
b- waves 


v  -  0.47  , 
u  -  0 . 11  , 
v  -  0 . 48 


(9.6) 


For  a  travelling  incident  wave  with  >  u,  there  are  two  resonances  as  a 


function  of  u,  and  for  an  evanescent  incident  wave  with  >  u  there  is  only  a 
single  resonance  as  a  function  of  v.  If  <  w  there  is  only  one  resonance  for 
a  travelling  incident  wave,  but  two  for  an  evanescent  wave.  If  p  *  1  and  7=0 
then  the  two  resonances  are  always  very  close  together.  This  can  be  understood 
from  the  fact  that  «  is  very  close  to  unity  in  this  limit. 

10.  Specular  reflection 

Ir  the  limit  of  a  weak  interaction  (7=0),  in  combination  with  close 

resonance  (p  =  1) ,  the  Fresnel  coefficients  for  the  specularly- reflected  waves 

at  the  incident  frequency  are  vanishingly  small  (section  8).  For  large 

2  2 

angles  of  incident  9,  the  parameter  u  -  cos  9  (travelling  incident  wave)  can 
be  on  the  order  of  or  |p-l|,  in  which  case  the  approximations  of  section  8 
are  not  accurate.  Since  the  specular  wave  is  also  reflected  back  into  the 
region  z  >  0,  it  will  interfere  with  the  pc-wave,  and  therefore  we  cannot 
neglect  this  component  in  the  situation  of  grazing  incidence.  Figures  7  and  8 
illustrate  the  behavior  of  |Xr|  for  various  angles  of  incidence.  For  7=0  and 
u  =  1,  |Xrj  indeed  disappears,  but  for  u  -»  0  the  value  of  |Xr|  approaches 
unity.  Phase  matching  between  the  a-  and  b-waves  in  the  PC  is  again 
responsible  for  the  oscillatory  behavior  of  |Xr|  as  a  function  of  d. 

11.  Conclusions 

We  have  studied  the  scattering  of  travelling  and  evanescent  waves  by  a 
phase  conjugator  in  the  four-wave  mixing  configuration,  without  restrictions  on 
the  angle  of  incidence,  the  interaction  strength  or  the  frequency  detuning  with 
the  pump  beams.  The  Fresnel  coefficients  for  the  various  waves  were  derived 
from  Maxwell's  equations,  subject  to  the  appropriate  boundary  conditions, 


without  the  usual  slowly -varying  amplitude  approximation.  It  was  shown  that  in 
the  limit  of  weak  coupling  and  perfect  resonance,  the  reflection  coefficient 
for  the  pc-wave  reduces  to  the  well-known  result  (8.8),  which  implies  the 
resonance  condition  (8.9)  for  a  four- wave  mixing  PC.  We  were  able  to  track 
down  the  origin  of  these  resonances  to  perfect  phase  matching  between  the  two 
(a  and  b)  modes  of  the  PC.  In  addition,  we  found  strong  resonances  at  these 
values  of  the  parameters  where  one  of  the  waves  (a,  b  or  pc)  turns  from  a 
travelling  wave  into  an  evanescent  wave.  We  were  not  able  to  find  a  convincing 
physical  explanation  for  these  resonances,  but  from  numerical  examples  it 
follows  that  they  are  definitely  present.  Finally,  we  showed  that  for  large 
angles  of  incidence  the  nonlinear  specularly-reflected  wave  has  an  amplitude  of 
the  same  order  as  the  pc-wave. 
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Figure  1.  Dispersion  relation  for  a  phase  conjugator.  Curves  a  and  b  are 

2  2  2  2  - 
ka/k  and  k^/k  ,  respectively,  as  a  function  of  and  the  coupling 

parameter  is  7^  -  0.2.  The  discontinuous  behavior  around  -  w  comes  from  the 

choice  of  the  roots  in  the  definition  of  n  and  n,  . 

'a  'b 

Figure  2.  Schematic  representation  of  the  various  waves  from  section  5.  The 
arrows  indicate  the  wave  vectors,  which  all  have  the  same  parallel  component 
ki| .  Their  perpendicular  components  are  approximately  the  same,  apart  from 
the  sign.  The  inc- ,  t-  are  r-waves  are  u^-waves  and  they  travel  in  the 
direction  of  the  wave  vector.  The  pc-  and  ni-waves  have  a  frequency  u> ^  <  0, 
and  therefore  they  travel  in  the  direction  opposite  to  the  wave  vector.  This 
is  indicated  by  a  circle  on  the  arrows.  Inside  the  PC  we  have  four  different 
wave  vectors  and  every  vector  corresponds  to  two  fields  according  to  Eqs.  (5.9) 
and  (5.10).  The  a+  and  a-  fields  are  essentially  u^-fields,  and  the  b+  and  b- 
fields  are  negative- frequency  waves.  For  7  0  these  principle  waves  couple  to 

a  field  with  the  same  wave  vector  but  with  a  frequency  of  opposite  sign.  These 
fields  are  indicated  by  broken  arrows.  The  wave  which  couples  to  the  principle 
wave  always  propagates  in  the  opposite  direction. 

Figure  3.  Absolute  value  of  the  reflection  coefficient  for  the  pc-wave  with  s- 
polarization  for  u  -  1  and  7  -  0.05,  as  a  function  of  the  layer  thickness  d. 

Curves  a,  b,  c  and  d  correspond  to  -  1.006,  1.01,  1.05  and  1.1, 

respectively.  The  peaks  result  from  interference  between  a-  and  b-waves,  and 
the  positions  of  the  peaks  are  given  by  Eq.  (8.11).  For  larger  detunings  the 

peaks  shift  and  become  less  pronounced,  in  agreement  with  Eq .  (8.12).  In  the 

limit  -  u>,  curve  a  turns  into  an  extremely  sharp  resonance. 
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Figure  4.  Reflectivity  for  the  phase -conjugated  wave  as  a  function  of  d.  The 
parameters  are  1q  -  0.05  and  w^/w  -  1.01.  Plot  a  corresponds  to  a  travelling 
incident  wave  (u  -  0.25),  and  curve  b  represents  an  evanescent  field  (u  - 
0.25i).  The  sharp  resonances  in  curve  a  originate  again  from  phase 
interference,  and  their  position  is  given  by  Eq.  (8.11).  It  is  seen  that  the 
reflectivity  for  an  evanescent  wave  is  almost  negligible,  as  compared  to  a 
travelling  wave. 


Figure  5 .  Absolute  value  of  X^c  as  a  function  of  u  for  7-0.05,  d  -  5  and 

-  1.05  (curve  a)  and  -  0.95  (curve  b) .  The  two  peaks  in  curve  a  are 

situated  at  ures  -  0.43  and  urgs  -  0.44,  and  they  appear  because  the  pc-wave 
and  the  b-wave  become  evanescent  for  lower  values  of  u.  In  curve  b  the  left¬ 
most  peak  corresponds  to  an  evanescent  a-wave  for  lower  values  of  u,  and  the 
other  peak  is  a  phase -matching  interference  from  section  8. 


Figure  6.  Same  as  Fig.  5,  but  now  as  a  function  of  v  -  -iu.  For  v  larger  than 
the  resonance  of  curve  a  the  a-wave  is  evanescent,  and  for  v  larger  than  0.48, 
both  the  pc-wave  and  the  b-wave  are  evanescent,  which  gives  rise  to  the  peaks 
in  curve  b. 


Figure  7 .  Reflection  coefficient  for  the  specular  wave  as  a  function  of  the 
normalized  layer  thickness  d.  The  parameters  are  w^/u>  -  1.01  and  7  -  0.05. 

For  curves  a,  b  and  c  we  took  u  -  0.2,  0.25  and  0.45,  respectively.  It  is  seen 
that  |Xr|  is  not  small  for  large  angles  of  incidence. 


Figure  8 .  Same  as  Fig.  7  but  with  u  -  0.01  and  0.1  for  curves  a  and  b, 
respectively.  In  the  limit  u  -»  0,  |Xr|  approaches  unity,  except  at  the  sharp 
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